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Abstract 

The e + e~ — > ZZ cross section at y/s = 182.7 and 188.6 GeV has been measured using 
the aleph detector. The analysis covers all of the visible ZZ final states and yields cross 
section measurements of 

or zz (182.7 GeV) = 0.11 ± (stat.) ± 0.04 (syst.) pb 

and 

(Tzz(188.6 GeV) = 0.67 ± 0.13 (stat.) ± 0.04 (syst.) pb, 
consistent with the Standard Model expectations. 
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1 Introduction 



The successful operation of lep at and above the ZZ threshold in 1997 and 1998 allows 
observation of a sizable number of pair-produced resonant Z bosons. Within the Standard 
Model, the process e + e~ — > TIL proceeds via the two MC2 diagrams which involve the 
t-channel exchange of an electron (see Figure [I]). A measurement of this process, the 
dominant irreducible background for Higgs boson searches at lep, tests the rates predicted 
by the Standard Model. 

This letter describes a measurement of the ZZ MC2 cross section at a centre-of-mass 
energy of 188.6 GeV using the aleph detector. It begins with a description of the detector 
and data samples, continues with a description of the event selections for the various final 
states, and concludes with the measured cross section. Application of the same analysis 
to the 1997 aleph data also yields a measurement of the cross section at y/s =182.7 GeV. 

2 ALEPH Detector 

This section briefly describes the most salient features of the aleph detector; detailed 
descriptions of the detector and its performance can be found elsewhere |2|]. 

aleph is a cylindrically symmetric detector with its axis coincident with the beam 
axis. The three innermost detectors — a silicon microstrip detector, the inner tracking 
chamber, and the time-projection chamber (TPC) — measure the momentum of charged 
particles. With a 1.5 T axial magnetic field provided by a superconducting solenoidal 
coil, these detectors together achieve a transverse momentum resolution of Sp t /pt = 
6 x 10" 4 ^ © 5 x 10~ 3 (p t in GeV/c). 

The excellent resolution of the silicon microstrip detector facilitates the identification 
of jets containing b hadrons. A neural network combines impact parameter and secondary 
vertex information with other jet information such as track multiplicity, track rapidity, 
and the presence of leptons, to provide a measure of the b content of jets. 

An electromagnetic calorimeter placed between the TPC and the superconducting 
coil identifies electrons and photons, and measures their energies with a resolution of 
5E/E = 0.18/-/E + 0.009 (E in GeV). This sampling calorimeter has a depth of 22 
radiation lengths and consists of projective towers each of which subtends a solid angle 
of approximately 0.9° x 0.9°. 

The iron return yoke is instrumented with 23 layers of streamer tubes and serves 
as a hadronic calorimeter. It is seven interaction lengths deep and achieves an energy 
resolution of SE/E = 0.85/ \/E for charged and neutral hadrons. Two additional layers 
of streamer tubes outside the return yoke aid the identification of muons. 

The luminosity is measured with small-angle Bhabha events, using lead-proportional 
wire sampling calorimeters ||, with an accepted Bhabha cross section of 4.6 nb at 
182.7 GeV and 4.3 nb at 188.6 GeV [|. Integrated luminosities of 56.81 ± 0.11 (stat.) ± 
0.29 (syst.) pb" 1 and 174.21 ± 0.20 (stat.) ± 0.82 (syst.) pb" 1 were recorded respectively, 
at mean centre-of-mass energies of 182.66 ± 0.05 GeV and 188.63 ± 0.04 GeV ||. 

An energy flow algorithm [[| combines the information from the tracking detectors 
and calorimeters and provides a list of reconstructed charged and neutral particles. These 
energy flow objects are used in the analysis described below. 
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Figure 1: The t-channel exchange diagrams show the e + e — > ZZ A/"Ci? processes in the 
Standard Model. 



3 Simulation of Signal and Background 

The yfszz Monte Carlo generator || provides a calculation of the expected Standard 
Model cross section for the J\fC2 processes. The expected cross section is 0.26 pb at 
y/s = 182.7 GeV and 0.65 pb at 188.6 GeV. With a consistent set of electroweak 
parameters, three other Monte Carlo generators — pythia J7|, excalibur [§], and 
GRACE4F |9| — all agree with the yfszz calculation within 2%. This value is taken as 
the theoretical uncertainty on the MC2 cross section determination. 

Samples of ZZ compatible events from the MC2 and MC8 diagrams were generated 
with the pythia Monte Carlo generator. The efficiencies were determined from the MC2 
sample. The MC8 diagrams include, in addition to the diagrams of Figure [l], six other 
diagrams in which one or both Z bosons are replaced by a virtual photon. The difference 
between the NC2 sample and the NC8 sample was taken as background. This procedure 
accounts for any interference between the four-fermion diagrams. As a check, the analysis 
was repeated with samples generated with the excalibur Monte Carlo generator (see 
Section |5TT| ) . 



The koralw ||T0[ Monte Carlo generator was used to create a sample of WW-like 
four-fermion background events. To avoid double counting, events in the MC8 TIL sample 
with a WW-like final state are removed. As the koralw generator does not include the 
part of the e + e~ — > Weu and e + e~ — > Ze + e~ processes with an electron near the beam 
axis and outside of the detector fiducial volume, additional samples for these processes 
were generated with the pythia Monte Carlo generator. 

For the two-fermion background, e + e~ — * and e + e~ — > t + t~ events were 



generated with koralz JTT|, e + e — > e + e events with bhwide [|12||, and e + e — > qq 



events with pythia. The two-photon background — 77 — > e + e~, 77 — > fi and 
77 — > t + t~ — was generated with PHOT02 JTB|. Studies showed that the process 
77 — > hadrons contributes negligibly to the background expectation after all of the 
analysis cuts. 
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Table 1: The efficiencies, ZZ signal expectation, background expectation, number of 
observed events, and the measured MC2 cross section with the statistical error for each 
channel at y/s = 188.6 GeV. 
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4 Event Selection 

The analyses described below select events from all visible ZZ final states. While based 



on the Higgs boson search [14], these analyses have been tailored to the ZZ production 
process. 

Two parallel analyses measure the ZZ production cross section. The first is an entirely 
cut-based analysis which uses the qqqq, qqz/17, ££XX, and llvv channels. Throughout this 
paper, the symbol t denotes an electron or muon and X denotes a quark or charged lepton. 
The second analysis, the "NN" analysis, replaces the cut-based selections in the qqqq and 
qqz/17 channels with neutral networks and introduces an additional neural network to 
identify the rrqq final state. 

Table [T] summarizes the results of the selections described below. 



4.1 Selection of ZZ — > UXX Final States 

This analysis identifies those events in which one of the Z bosons decays into a pair of 
electrons or muons and the other decays into anything except neutrinos. Although the 
branching fraction of Z — > tl is small, the analysis benefits greatly from the high lepton 
identification efficiency and excellent mass resolution. 

Selected events have four or more reconstructed charged particles. The total energy 
of all the charged particles must exceed 10% of the centre-of-mass energy. To remove 
radiative returns to the Z, the most energetic isolated photon in an event must have an 
energy less than 0.75 7 pea k, where 7 pea k is the most likely energy of the radiative return 
photon and is given by 7 pea k = 0.5 [yfs — m|/ \fs). 



Standard lepton cuts are used to identify electrons and muons in the event [ D| . 
Electron energies are corrected for possible bremsstrahlung losses. To maintain a high 
selection efficiency, isolated charged particles are also considered as lepton candidates. 
The isolation is defined as the half-angle of the largest cone about a particle's direction 
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which contains 5% or less of the event's total energy. Isolated particles have an isolation 
angle greater than 10°. The pair of oppositely-charged lepton candidates with an invariant 
mass closest to the Z mass is chosen as the Z — > ii pair. Only those pairs with at least 
one identified lepton are considered and efi combinations are rejected. Photons consistent 
with final state radiation are included in the calculation of the invariant mass of the 
lepton pair. Events in which the lepton pair contains an electron consistent with one 
from a photon conversion are rejected. 



After the Z — > ££ lepton pair has been selected, the DURHAM algorithm |TJ| clusters 
the remainder of the event into two jets; these jets must contain at least one charged 
particle. To remove Z7* events, the selection requires that the invariant mass of the two 
jets exceed 15 GeV/c 2 . To remove qq events, the sum of the transverse momenta of the 
leptons with respect to the nearest jet must be greater than 20 GeV/c, where the nearest 
jet is the one which forms the smallest invariant mass with the lepton. 

The process WW — > qjqiv constitutes a large background for the selected events in 
which the Z — > 11 pair contains one identified lepton and one isolated charged particle. 
The two W masses are calculated assuming that the event comes from this process. The 
mass mi ep is the invariant mass of the identified lepton and the missing momentum; the 
other mass m^d is the invariant mass of the remainder of the event. For selected events, 
these masses satisfy either m had + mi ep < 150 GeV/c 2 or m had — mi op > 20 GeV/c 2 . 

The mass of the lepton pair is used as the mass of the first Z candidate; that of the 
second candidate is the mass recoiling against the two leptons. Requiring the above cuts 
and that both Z masses be larger than 30 GeV/c 2 , the analysis selects 92 events from the 
data with 90.6 ± 1.4 expected from signal plus background. Figure |2| shows the dilepton 
invariant mass distribution. 

To further reduce background levels, the two reconstructed masses must be consistent 
with the Z boson mass. As the leptonic invariant mass and recoil mass have different 
resolutions, an elliptical cut is defined using 



.2 



' m u - m z \ 2 ^ ( m recoi i - m z \ 2 

°m« / \ °m rccoil / 



where a mu = 2.5 GeV/c 2 and cr mrccoil = 3.3 GeV/c 2 . Selected events have r < 3. An ££XX 
candidate, classified as a ZZ — > fifirr event, is shown in Figure ^. 

Systematic uncertainties from lepton identification, tracking resolutions, and 
reconstruction of kinematic variables have been studied. The total relative systematic 
uncertainty on the efficiency is 1.3% and the relative uncertainty on the expected number 
of background events is 27%. Both uncertainties are dominated by limited Monte Carlo 
statistics. 



4.2 Selection of ZZ — ► qqqq Final States 

A total reconstructed energy near -y/i and four distinct jets characterize events from the 
ZZ — > qqqq final state. This final state has the largest branching fraction, but also the 
largest background. 

The cut-based and NN-based analyses of this channel use a common preselection. The 
event must contain eight or more charged particles and the total energy of all charged 
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Figure 2: The distribution of the dilepton masses for the MXX channel. All cuts except the 
elliptical mass cut have been applied. The grey histogram shows the expected background 
and the hollow histogram shows the expected signal. 

particles must exceed 10% of the centre-of-mass energy. 

The DURHAM algorithm clusters the event into four jets. To remove those events 
inconsistent with a four-jet topology, the value of y 34: , the y cut in the clustering algorithm 
which changes the event from a three-jet to four-jet topology, must be larger than 0.004. 
Additionally, each jet must contain at least one charged track. 

Two additional cuts suppress events with initial state radiation. First, the longitudinal 
momentum of the event must satisfy the requirement \p z \ < 1.5 (M vis — 90) with M vis in 
GeV/c 2 and p z in GeV/c, to remove events in which the initial state photon is along 
the beam axis and outside of the detector fiducial volume. Second, the electromagnetic 
fraction of all of the jets must be smaller than 80% to remove those events in which the 
photon is inside the fiducial volume. This second cut also removes the qqe + e~ overlap 
with the leptonic channel described above. 

4.2.1 Cut-based Selection 

The cut analysis tightens the preselection with three additional cuts. The first two target 
qq events: the thrust value of the event must be smaller than 0.9 and the sum of the four 
smallest angles between jets must be greater than 350°. The third cut reduces the overlap 
of qq/i + /i~ events between the qqqq and the ££XX selections. Accepted events have an 
invariant mass of the two most energetic muon candidates less than 50 GeV/c 2 and 
satisfy the condition pi+p? — "m^c < 35 GeV/c, where p± and P2 are the momenta of the 
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Figure 3: An event selected by the ££XX analysis and classified as a ZZ — ► /j,/! tt candidate. 
Tracks 2 and 4 form the Z — > candidate with an invariant mass of 96.0 GeV/c 2 . Tracks 
1 and 3 form the Z —>■ tt candidate with a recoil mass of 90.3 GeV/c 2 . 



two muons. For events with only one identified muon, this cut reduces to p± < 35 GeV/ c. 
These negligibly affect the efficiency for the qqqq final state but reject more than 97% of 
the qq/i + /i~ overlap and also a small fraction of the WW background. 

After the above preselection, 1352 ± 6 events are expected from the simulation, 
approximately 80% of which are WW — > qqqq events. In the data, 1219 events are 
observed. This discrepancy has already been reported by aleph in a preliminary 
measurement of the WW — > qqqq cross section [[17]]. Its impact on the ZZ cross section 
measurement is studied in Section |6~2" . 

The selection subsequently relies on the tagging of jets from b quarks and on the dijet 
mass information. A four-constraint fit [18| is then applied to improve the jet momentum 
resolution. A b-tagging neural network produces a value rji for each jet which is near 
unity for b jets and near zero for other jets fl4l . 

The four-jet analysis is split into three sub-analyses: a bbbb selection, a bbqq selection, 
and a non-b quark selection. When calculating the cross section, the two b selections are 
treated as a single channel. The last sub-analysis replaces the b-tagging cuts with strict 
mass requirements to reject background while retaining efficiency for the majority of ZZ 
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events (~62%) which do not contain b quark jets. 

The bbbb analysis selects events with high b content, well-isolated jets, and large dijet 
masses. It requires that yu > 0.020, that the sum of the dijet masses for at least one of 
the dijet combinations be above 170 GeV/c 2 , and that 9.5 yu + J2Vi > 3.1 where the sum 
is over the four reconstructed jets. 

For the other sub-selections the b-tagging requirements give way to tighter dijet mass 
requirements. The mass information is contained in the quantities xw an d Xz defined as 

2 (m 12 + m 34 - 2mA 2 ( m 12 - m u \ 2 

where i stands for W or Z. The quantities a l s and a l D are the resolutions of the sum and the 
difference of the dijet masses for the correct dijet combination, respectively. The contours 
of constant xw and xz define two ellipses referred to as WW and ZZ ellipses throughout 
this section. For the WW and ZZ ellipses respectively, the values 4 and 3 GeV/c 2 are used 
for a l s and the values 10 and 16 GeV/c 2 , for a 2 D . 

The bbqq selection requires that at least one dijet combination falls inside the ZZ 
ellipse with xz < 2.40. For that combination, the dijet not containing the most poorly 
b-tagged jet must be compatible with Z — > bb in terms of b-tagging: min(^i,^ 2 ) > 0.20 
and - log 10 (l - r/i)(l - r] 2 ) > 1.50. 

Finally, the selection for qqqq events without b jets raises the 2/34 cut to 0.006, tightens 
further the ZZ ellipse to xz < 1-75, and reduces the WW background by requiring that no 
dijet combination falls inside a WW ellipse with xw < 1-60. To maintain the statistical 
independence of the b and non-b selections, this selection accepts only those events which 
have not been previously selected by the bbbb or bbqq cuts. 

Figure |] shows the sum of the dijet masses for the bbbb and bbqq selections with the 
elliptical mass cut removed. 

Although this selection does not specifically target rrqq events, it nonetheless has a 
significant efficiency for these events (~17%). Consequently, the rrqq selection described 
below is not included in the cut-based cross section measurement. 

4.2.2 Neural Network Selection 

A multivariate neural network is trained to select ZZ signal events from qq and WW 
background events. During training, the network is presented with up to six dijet 
combinations per qq or WW event, but only the correct dijet pairing for each signal event. 
If one Z boson decays to b quarks, those jets are labelled 3 and 4, and the other two jets 
are labelled 1 and 2. Otherwise the jet labels are random. Reconstructed dijet masses 
are obtained using a four-constraint fit which requires energy-momentum conservation. 
Twenty-three event variables are used to discriminate specific signal and background 
features. In many cases, inclusion of similar variables in different combinations improves 
the neural network training. 

The event variables 2/34, thrust, and sphericity aid in selecting events with four 
separated jets. The missing energy and max (i? charged track /_Ej Ct ), the maximum scaled 
energy among all jets of the most energetic charged track in a jet, specifically reject 
semileptonic W decays. 
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Figure 4: The distribution of the sum of the dijet masses (m 12 + m u ) for the bbbb 
and bbqq channels. All cuts have been applied except the elliptical mass cut. The hollow 
histogram is the contribution from the signal and the darker region that from the expected 
background. 

A b-tagging neural network |L4| gives an output rji for each jet, and these variables 
are input as min^, 774), (1 - ry 3 )(l - 774), and Y$=i Vi- 

The reconstructed dijet masses aid the selection of resonant ZZ production events. The 
mass combinations used as network inputs are m 12 , m 34 , and (mi 2 — m^) 2 + (m 3i — mz) 2 . 
These variables perform the same function as an elliptical mass cut. 

The jet boosted sphericity (calculated in the rest frame of the jet) of the two jets with 
lowest rji, the track multiplicity in these jets, and the two lowest jet masses all help to 
separate light quark jets from gluon jets. The track multiplicity in this case is restricted 
to tracks with a rapidity with respect to the jet axis greater than 1.6. 

The sum of the four smallest inter-jet angles discriminates against qq background 
events. The angular variable min (cos^ + cos 6*^) minimized over all jet combinations 
selects events with pairs of back-to-back jets. The event broadening variable B = 
min(B + , B_) also aids in this discrimination by rejecting qq background with gluons 
that can lead to broad jets. The broadening variable for each of the two hemispheres S± 
is given by B± = (J2ieS±PTi)/(J2i \Pi\), where is calculated with respect to the thrust 
axis and the sums run over tracks. Finally, the largest jet energy and the two smallest jet 
energies improve the overall discriminating power. 

Of all of the dijet combinations, the one with the largest neural network output is 
used. The output of the neural network is shown in Figure |5] for Monte Carlo simulation 
and for data. Only events with a network output greater than 0.7 are used in the cross 
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Figure 5: Neural network output from the qqqq selection for data and Monte Carlo 
simulation. The grey histogram shows the expected background, and the hollow histogram 
shows the expected MC2 signal. 



section calculation. 



4.2.3 Systematic Uncertainties 

The systematic uncertainties in this channel include effects from the modelling of 
b physics, from discrepancies in tracking between the simulation and the data, from 
discrepancies in reconstructed jet kinematics, from uncertainties in the background cross 
sections, and from uncertainties in gluon splitting into heavy flavours. 

For the cut-based analysis, the total relative systematic uncertainties on the signal 
efficiencies are 1.7% and 1.3% for the b and non-b selections, respectively, with all of the 
sources having comparable contributions. The relative uncertainties for the background 
are 13% and 4%, dominated by the limited Monte Carlo statistics and jet corrections. 

For the NN-based analysis, the total relative systematic uncertainty on the signal 
efficiency is 3%. No single effect dominates the uncertainty. The total relative 
uncertainties for the background are 11% (qq) and 7% (WW), dominated by uncertainties 
in the reconstruction of kinematic variables and limited Monte Carlo statistics. 



4.3 Selection of ZZ — > qqyi/ Final States 

Approximately 30% of the ZZ events have a qqz/z7 final state characterized by a missing 
mass $f and visible mass consistent with the Z mass. 
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Both the cut-based and NN-based analyses for this channel share a common 
preselection. Events must contain more than four reconstructed charged particles and the 
total energy of all the charged particles must exceed 10% of the centre-of-mass energy. 
The plane perpendicular to the thrust axis divides the event into two hemispheres. These 
hemispheres are the "jets" used to calculate all of the kinematic quantities below. Both 
of the hemispheres must have a nonzero energy. 

The preselection rejects events with a total energy more than 30° degrees from the 
beam axis smaller than 0.25 t/s and a missing transverse momentum smaller than 0.05 y/s. 
This cut primarily removes 77 — > hadrons events. 

Two additional cuts remove much of the background from qq events with initial state 
radiation. First, the magnitude of the longitudinal event momentum must be less than 
50 GeV/c. Secondly, the missing mass must be greater than 50 GeV/c 2 . 

After the preselection, the background is dominated by qq and WW events. Assuming 
the Standard Model cross section for ZZ production, 1533 ±7 events are expected in good 
agreement with the 1493 actually observed in the data. 



4.3.1 Cut-based Selection 

The two-fermion background remaining after the preselection consists largely of qq events 
with two or more initial state radiative photons. Requiring the acoplanarity A 12 to be 
greater than 0.08 removes much of this background. The acoplanarity of the two jets is 
A12 = (ji x J 2) ' £ where % are the unit vectors along the jet directions. 

Pairs of W bosons which subsequently decay as WW — > tv qq form the other major 
background. To reject cases where the tau lepton decays leptonically, a cut on the isolation 
of identified leptons is made. The isolation is defined as the sum of the energy in a 30° 
cone around the lepton direction. Selected events have an isolation greater than 13 GeV. 
When the lepton from the tau decay is unidentified or the tau lepton decays hadronically, 
the previous cut is ineffective. To remove these cases, the analysis requires that the angle 
between the reconstructed charged particle with the highest momentum and the nearest 
charged particle be less than 20°. Additionally, the event is reclustered into minijets with 
the jade algorithm [[TjJ using y cut = (2 GeV/y/s) 2 . The energy of the most isolated 
minijet must be smaller than 8 GeV. 

To remove events from the e + e~ — > Wez/ and e + e~ — > Ze + e~ processes which have 
a detected final state electron near the beam axis, the analysis requires that the energy 
in a cone of 12° around the beam axis must be smaller than 0.02 y/s. To remove events 
which have a missing momentum near the beam axis and may be poorly measured, the 
smallest angle between the missing momentum and the beam axis must be larger than 
25°. Figure ^| shows the distribution of the reconstructed mass for candidates which pass 
all of the above cuts. The reconstructed mass m rec is the invariant mass of the two jets 
calculated with the missing mass constrained to mz- 

Because the mass of the Z boson is known, the mass information can be used effectively 
to reject much of the background. As for the other channels, this selection uses an elliptical 
mass cut of "radius" 

r 2 _ ( m rcc -m z \ 2 ( W -m z \ 2 
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Figure 6: The distribution of reconstructed masses for the qqz/z7 channel. The grey region 
is the expected contribution from the background and the hollow region that from the ZZ 
MC2 signal. 

where cr mrec = 3.1 GeV/c 2 and = 8.5 GeV/c 2 are the resolutions. Selected events have 
r < 2. 

4.3.2 Neural Network Selection 

After preselection, a 12-variable neural network is employed to distinguish signal events 
from background. Four of the variables used in the neural network analysis (acoplanarity, 
minijet energy, direction of missing momentum, and energy within a 12° cone around the 
beam axis) are shared with the cut analysis. Additionally, the neural network analysis 
uses the reconstructed Z boson mass and the missing mass directly rather than using an 
elliptical mass cut. This analysis also includes two b-tagging variables although only a 
small fraction of events have a bbz/77 final state. In addition, acollinearity, the visible 
energy beyond 30° of the beam axis, the energy within a 30° wedge of the direction of the 
missing momentum, and the z-component of the momentum are used as variables. The 
output of the neural network is shown in Figure [7] for Monte Carlo simulation and for 
data. Only events with a network output greater than 0.5 are used for the cross section 
calculation. 

4.3.3 Systematic Uncertainties 

The results have been corrected for unsimulated accelerator background events which 
increase the energy deposited near the beam axis; this effect was studied using random 
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Figure 7: Neural network output from the qqW selection for data and Monte Carlo 
simulation. The grey histogram shows the expected background, and the hollow histogram 
shows the expected MC2 signal. 



triggers. In addition, half of the correction is used as an estimate of the systematic 
uncertainty. Similarly, the jet kinematics are corrected to make them correspond better 
to those observed in the data. Again, half of this correction is used as the systematic 
uncertainty. The systematic uncertainties also include uncertainties from the background 
production cross sections and the limited statistics of the Monte Carlo samples. 

Adding these uncertainties in quadrature gives total relative uncertainties of 2% and 
5% on the signal and background, respectively, for the cut-based analysis. For the 
efficiency, all of the sources are significant; for the background, the limited Monte Carlo 
statistics dominate. 

For the NN-based analysis, the b-tagging systematic uncertainties are estimated by 
reweighting the data with a function obtained from comparison of Z-peak Monte Carlo 
simulation and data. The relative systematic errors on signal and background are 2% and 
3%, respectively. 



4.4 Selection of ZZ — > llvu Final States 

The preselection requires events to have exactly two identified electrons or muons and 
no other reconstructed charged particles. These leptons must have the same flavour and 
opposite charge. The fraction of the centre-of-mass energy deposited in the detector more 
than 30° from the beam axis ^30° must satisfy the requirement 0.4 < / 30 o < 0.6. The 
angle between the two lepton directions must be less than 178°. For electron pair events, 
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49 ± 2 events are expected and 43 observed; for muon pair events, 13.2 ± 0.3 events are 
expected and 10 observed. 

Cuts on other kinematic quantities further suppress the background. The invariant 
mass of the leptons and the missing mass are required to be within an ellipse of r < 1.7 
defined as above with a pi = 3.3 GeV/c 2 and a mu = 2.5 GeV/c 2 . The angle between the 
missing momentum and the beam axis must be larger than 6.7°, and the total energy not 
associated with the leptons must be less than 5.6 GeV. 

The systematic uncertainty is 5% for the efficiency and 17% for the background. These 
include uncertainties from lepton identification and tracking resolution but are dominated 
by the limited Monte Carlo statistics. 

4.5 Selection of ZZ — > rrqq Final States 

In the neural network based selection, a dedicated analysis selects ZZ — > rrqq events with 
a higher efficiency than the cut-based four-jet selection. 

Events must contain at least eight tracks and the total energy of all the charged 
particles must exceed 0.20 ^/s. Background from WW and ZZ events with electron or 
muon decays is suppressed by rejecting events having an identified lepton with energy 
greater than 0.25 y^s. Radiative returns to the Z peak, characterized by high missing 
energy and high missing longitudinal momentum, are rejected by requiring \p z \+ ft < 
1-8 7 P eak, where p z and $ are respectively the missing longitudinal momentum and missing 
energy. To further reject radiative returns, events are also rejected if the event \p z \ is 
greater than 0.6 7 pca k- 

Events passing the preselection cuts are clustered into minijets having invariant masses 
consistent with a tau lepton, m- iCt < 2.7 GeV/c 2 . The r candidates are selected from these 
minijets using a series of quality cuts based on multiplicity, isolation, and momentum. 
To be considered a r candidate, a minijet must have one, two, or three charged tracks 
with momenta larger than 1 GeV/c. If the minijet has three charged tracks, it must be of 
unit charge; if the minijet has two charged tracks, the minijet charge is the charge of the 
track with higher momentum. The minijet isolation angle must be larger than 15°. This 
isolation angle is the half-angle of the largest cone around the minijet direction containing, 
in addition to the minijet energy, no more than 5% of the total event energy outside the 
cone. Finally, the energy of a two- or three-prong minijet must be greater than 12.5 GeV, 
while a one-prong minijet with less than 80% of its energy from charged particles must 
have an energy greater than 7.5 GeV. Otherwise, no momentum cut is made. 

Only events with at least two r candidates are further considered. At least one of the 
r minijets must have exactly one prong, and the two minijets must have opposite charge. 
The rest of the event is clustered into two jets using the Durham algorithm. All four jets 
in the event are rescaled using an overconstrained kinematic consistency fit in which the 
jet directions are fixed and the minijet masses are set to m T . The fit x 2 is calculated from 
energy-momentum conservation, the hadronic jet resolutions, and the difference between 
the two fitted dijet masses. This has the effect of constraining the dijets to equal masses. 
In no case are the hadronic jets allowed to rescale to less than 75% of their measured 
momenta. A typical event may have several possible combinations of potential r minijet 
candidates; only the combination with the smallest calculated kinematic \ 2 is further 
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considered. 

A five-variable neural network selects ZZ events from preselected events. Because 
the r leptons decay with one or two neutrinos, the event transverse momentum pt is 
input to discriminate against background events without missing energy. The \ 2 provides 
a measure of the event's kinematic consistency with the ZZ signal. The sum of the 
t candidates' isolation angles and the sum of the r candidates' transverse momenta 
with respect to their nearest hadronic jet distinguish r minijets, which are expected 
to be well-isolated, from fake r minijets on the periphery of a hadronic jet. Finally, 
the reconstructed Z mass, when combined with the implicit equal mass constraint, helps 
discriminate between MC2 events and background. Events with a network output greater 
than 0.77 are selected for the cross section calculation. 

The systematic uncertainties arise from uncertainties on the reconstructed jet energies 
and angles, the background cross sections, and the limited Monte Carlo statistics. The 
relative systematic errors for signal and background are 3.0% and 10.8%. The error for 
signal is distributed equally among the sources, while the background error is dominated 
by limited Monte Carlo statistics. 

5 Combination of Channels 

Table [l] summarizes the efficiencies, numbers of expected signal and background events, 
the number of observed events, and the measured cross section for each channel. A 
maximum likelihood fit determines the MC2 cross section for the e + e~ — * ZZ process by 
combining the information from all of the channels. 

For the cut-based analysis, only the total numbers of events are used in the likelihood 
combination. Overlaps between the various channels are less than 0.2% and are negligible. 

For the given signal and background numbers, toy simulated experiments are 
generated, and the expected uncertainty on the cross section measurement is derived 
from the width of the extracted cross section distribution. The expected relative statistical 
uncertainty on the 188.6 GeV cross section measurement is 19.7%. 

For the neural-network based analysis, a binned likelihood fit is used for the qqqq 
and qqVF channels where the shape of the neural-network distribution is used; for the 
other channels only the total numbers of events are used in the likelihood fit. The qqqq 
distribution is divided into six bins between 0.7 and 1.0, while the qqz/z7 distribution is 
divided into ten bins between 0.5 and 1.0. 

For the neural-network based analysis, the overlaps between channels are explicitly 
removed. About 0.2% of the events in the qqqq selection are rrqq events in which 
a three-prong r is misidentified as a quark jet, and about 0.5% are qqu + /i~ events. To 
ensure the exclusivity of the analyses, the overlapping signal and background are explicitly 
subtracted from the qqqq results. The expected uncertainty after combining all channels 
is 18.0%. 

The systematic uncertainties were determined by adding a Gaussian smearing to the 
efficiencies and background estimates in toy Monte Carlo experiments and observing the 
change in the total error. Correlated contributions to the total systematic uncertainty of 
the measurement include the luminosity (0.5%) and the uncertainties in the WW and qq 
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Figure 8: The measured MC2 cross section as a function of the centre-of-mass energy 
compared to the expected Standard Model cross section. The shaded area represents the 
theoretical uncertainty (2%) on the yfszz calculations. 

cross sections. The total relative systematic uncertainties are 4% and 6% for the cut-based 
and neural-network based analyses, respectively. 

The measured cross sections at y/s = 188.6 GeV are 

a zz = 0.69 ± ° il (stat.) ± 0.03 (syst.) pb 

and 

azz = 0.64 ± ° il (stat.) ± 0.04 (syst.) pb 

for the cut-based and neural-network based analyses, respectively, compared to the 
Standard Model expectation of 0.65 pb. Table p] shows the cross section for each channel 
individually. 

The cut-based analysis selects 13 candidates with an expected background of 12.2 
events when applied to the 1997 data sample. The total efficiency for all channels is 33%. 
Combination of the channels with the above likelihood function yields a measured cross 
section of 

a zz = 0.11 ± ° il (stat.) ± 0.04 (syst.) pb 

at a centre-of-mass energy of \fs = 182.7 GeV, compared to the Standard Model 
expectation of 0.26 pb. The measurement is also consistent with a similar measurement 



from the L3 collaboration 20 . 
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6 Cross Checks 



6.1 Four-fermion Interference 

The pythia Monte Carlo generator does not fully simulate the interference between 
all of the four-fermion diagrams. This potentially biases the cross section measurement 
presented here. Using excalibur, which does include the interference, to generate 
samples for the signal and four-fermion background and rerunning the cut-based analysis 
yields a relative difference in the cross section measurement of —1%. This is a negligible 
effect when compared to the magnitude of the statistical error. 

6.2 Correlations Between the WW and ZZ Cross Sections 

The preliminary measurement of the WW cross section in the qqqq channel performed 
by aleph on the 1998 data differs from the Standard Model by —10%. Although this 
could be a statistical effect, the effect of a real reduction of the WW cross section on the 
extracted ZZ cross section was investigated. As a check, the expected WW background of 
the four-quark channels of the cut-based analysis was reduced by 10% and the channels 
recombined. This leads to a +3.5% relative change in the measured ZZ cross section. 

7 Conclusions 

Two analyses have measured the ZZ cross section using the data taken at ^Js = 188.6 GeV 
in 1998. The results agree and are consistent with the Standard Model expectation of 
0.65 pb. The ZZ production cross section is taken to be the arithmetic average of the 
above two measurements 

(x zz (188.6 GeV) = 0.67 ± 0.13 (stat.) ± 0.04 (syst.) pb . 

Additionally, the cut-based analysis has been applied to the 1997 data sample yielding a 
measurement of 

ctzz(182.7 GeV) = 0.11 ± g;^ (stat.) ± 0.04 (syst.) pb . 
The measured cross section is compared to the Standard Model expectation in Figure |8|. 
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